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ABSTRACT 

The low-luminosity active galactic nuclei in NGC 3147, NGC 4203, and NGC 4579 have been imaged at four 
frequencies with the Very Long Baseline Array. The galaxies are unresolved at all frequencies, with size upper 
limits of lO^'-lO'* times the Schwarzschild radii of their central massive black holes. The spectral indices between 
1.7 and 5.0 GHz range from 0.2 to 0.4; one and possibly two of the galaxies show spectral turnovers between 5.0 
and 8.4 GHz. The high brightness temperatures (7i,>10^ K) and relatively straight spectra imply that free-free 
emission and/or absorption cannot account for the slightly inverted spectra. Although the radio properties of the 
cores superficially resemble predictions for advection-dominated accretion flows, the radio luminosities are too 
high compared to the X-ray luminosities. We suggest that the bulk of the radio emission is generated by a compact 
radio jet, which may coexist with a low radiative efficiency accretion flow. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: individual (NGC 3147, NGC 4203, 
NGC 4579) — galaxies: nuclei — radio continuum: galaxies 



1. INTRODUCTION 

Most nearby normal and active galaxies with bulges appear 
to have supermassive black holes (BHs) at their centers (see Ko- 
rmendy & Gebhardt 2001 for a review), although galaxies with- 
out bulges may lack BHs (Gebhardt et al. 2001). The central 
BHs are thought to power the low-luminosity active galactic nu- 
clei (LLAGNs) in local Seyfert galaxies and LINERs. In a num- 
ber of these LLAGNs, the radio emission is compact on mil- 
liarcsecond (mas) scales (Falcke et al. 2000; Nagar, Wilson, & 
■ Falcke 2001) when imaged with the Very Long Baseline Array 
(VLBA; Napier et al. 1993)^. They often have mildly inverted 
centimeter radio spectra with spectral indices a « +0.3, where 

oc v"" (Wrobel & Heeschen 1991; Slee et al. 1994; Nagar et 
al. 2000). In our recent VLA study of the Palomar Seyfert sam- 
ple (Ho, Filippenko, & Sargent 1997), we found that about half 
of the weak Seyferts selected by optical spectroscopy have flat 
or inverted radio spectra (Ho & Ulvestad 200 1 , hereafter HUO 1 ; 
Ulvestad & Ho 2001), These galaxies differ from more lumi- 
nous Seyferts, which typically have optically thin synchrotron 
spectra with a w -0.7 (Ulvestad & Wilson 1984, 1989; Ulves- 
tad & Ho 2001). Possible explanations discussed by Ulves- 
tad & Ho (2001) included free-free emission/absorption or an 
advection-dominated accretion flow (ADAF; see Narayan, Ma- 
hadevan, & Quataert 1998; Quataert 2001). 

The Galactic Center source Sgr A* also has a mildly inverted 
spectrum at centimeter wavelengths, with a « +0.3 from 1 GHz 
to over 100 GHz (Falcke et al. 1998). This spectrum is some- 
times attributed to an ADAF (e.g., Narayan, Yi, & Mahadevan 
1995; Manmoto, Mineshige, & Kusunose 1997). In ADAFs, 
synchrotron radiation from hot electrons in a > 10^ K gas is 
self-absorbed, yielding a radio spectral index between +0.3 and 
+0.4 (Mahadevan 1997; Narayan et al. 1998). X-ray emission 
arises from Compton scattering or bremsstrahlung, but is far 
below the Eddington limit, since most of the plasma's energy is 
advected into the BH by the ions. The transition from a classical 



thin accretion disk to an ADAF takes place within ^ 10^-10"* 
Schwarzschild radii (/?s) of the BH. Compact radio emission 
from a number of LLAGNs has been attributed to ADAF pro- 
cesses similar to those inferred in Sgr A* (e.g., Fabian & Rees 
1995; Mahadevan 1997; Ho et al. 2000; Wrobel & Herrnstein 
2000). However, Di Matteo et al. (1999) find that the simplest 
ADAF models overpredict the radio and submillimeter fluxes 
in several nearby ellipticals. 

Another possible model for Sgr A* and LLAGNs is that of a 
compact radio jet in combination with an ADAF (Falcke 1996; 
Falcke & Markoff 2000; Yuan 2000). As reviewed by Falcke 
(2001), the jet accelerates slightly due to a longitudinal pres- 
sure gradient, and integration of jet radio emission may result 
in a slightly inverted spectrum. Compton upscattering of the 
radio photons produces the X-rays. An energetically significant 
jet or outflow will reduce the accretion rate in the innermost 
region of the ADAF, which may explain the low level of ra- 
dio emission seen in some ellipticals (Di Matteo et al. 1999; 
Quataert & Narayan 1999). 

In this Letter, we report four-frequency VLBA imaging of 
three LLAGNs imaged previously at arcsecond resolution with 
the VLA (HUOl). These galaxies were selected because they 
have inverted VLA spectra between 1.4 and 4.9 GHz and flux 
densities sufficiently high (>10 mJy) for easy imaging with the 
VLBA. Therefore, they were good candidates for testing the 
various models that might account for their spectral shapes. 

2. OBSERVATIONS AND RESULTS 

The VLBA observed NGC 3147, NGC 4203, and NGC 4579 
in March 2001 at central frequencies of 1.667, 2.271, 4.995, 
and 8.421 GHz. About six hours was spent integrating on each 
galaxy, divided among the four frequencies. Amplitude calibra- 
tion used a-priori gain values together with system temperatures 
measured during the observations, and typically is accurate to 
5%. Initial clock and atmospheric (phase) errors were derived 
from calibrator sources within 2°-A° of the galaxies. 
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TABLE 1 
VLBA Observational Results 



Galaxy 
(1) 


UT Date 
(2) 


a(J2000) 
{h m s) 
(3) 


<5(J2000) 

C ' ") 

(4) 


5l.7 

(mJy) 
(5) 


52.3 

(mJy) 

(6) 


55.0 

(mJy) 
(7) 


58.4 

(mJy) 
(8) 


Tb 

(K) 

(9) 


Calibrator 
(10) 


NGC 3147 


2001MarlO 


10 16 53.6511 


+73 24 02.695 


7.0 


6.8 


8.8 


9.3 


> 8.0 X 10** 


J1031+7441 


NGC 4203 


2001Marl9 


12 15 05.0553 


+33 11 50.384 


6.2 


7.4 


9.5 


10.0 


> 7.0 X 10** 


J1215+3448 


NGC 4579 


2001Mar26 


12 37 43.5222 


+11 49 05.488 


18.3 


21.1 


22.8 


19.3 


> 8.4 X 10** 


J1230+1223 



Note.— Col. (1) Galaxy name. Col. (2) UT observing date. Col. (3)-(4) VLBA position. Col. (5)-(8) Flux densities at 
1.7, 2.3, 5.0, and 8.4 GHz, respectively. Col. (9) Brightness temperature of the core. Col. (10) Calibrator source. 

NGC 3147 8.4 GHz NGC 4203 S.4 GHz NGC 4579 8.4 GHz 




MAS MAS 



Fig. 1.— 8.4 GHz VLBA images of (a) NGC 3147, (b) NGC 4203, and (c) NGC 4579. AH images are 9.6 mas on a side, with contour intervals increasing by 
factors of \/2 from 0.25 mJy beam"^ (negative contours are shown dashed). The restoring beam for each image is shown in the lower-right comer. 



using the technique of phase-referencing (Beasley & Conway 
1995). This resulted in galaxy position errors of ^ 1 mas, dom- 
inated by uncertainties in the calibrator positions (from John- 
ston et al. 1995; A. J. Beasley et al., in preparation). The initial 
phase calibration was limited at low frequencies by the active 
ionosphere, and at high frequencies by the wet troposphere, but 
each galaxy was detected at all four frequencies. Phase self- 
calibration was applied iteratively, resulting in noise-limited 
images. Approximate beam widths ranged from ~ 7 mas at 
1.7 GHz to 1.2 mas at 8.4 GHz, while the rms image noises 
were 90-150 /iJy beam"'. 

Gaussian fitting showed that all VLBA images were domi- 
nated by a single unresolved component; we estimate size up- 
per hmits of half the beam widths. Derived flux densities are 
accurate to ~ 5% at 2.3 and 5.0 GHz, and ~ 10% (due to larger 
self-calibration corrections) at 1.7 and 8.4 GHz. Table 1 sum- 
marizes the observations and main results, and Table 2 lists 
some derived quantities and other pertinent information. Fig- 
ure 1 shows the high-resolution, unresolved 8.4-GHz images of 
the galaxies, while Figure 2 compares the radio spectra with 
thatof Sgr A*. 

3. INDIVIDUAL GALAXIES 

NGC 3147.— NGC 3147 contains a Seyfert 2 nucleus in an 
She galaxy (Ho et al. 1997). The hard X-ray (2-10 keV) lu- 
minosity measured with A5CA is 3.4 x lO'^'erg s"' (Terashima, 
Ho, & Rak 2000a), and observations with the ROSAT HRI 

Gebhardt et al. (2000) use a = ae, the projected, luminosity-weighted velocity dispersion measured within the effective radius of the bulge. Since our objects do not 
have measurements of ae, we use a = ag, the central velocity dispersion; Gebhardt et al. (2000) show that ae ^ ag within a scatter of ~10%. 



(Roberts & Warwick 2000) show that this source is almost cer- 
tainly located at the nucleus. Our upper hmit to the source 
radius is -0.05 pc, or - 1500/?s for a BH mass of 3.6 x 10*^ 
Mq. We estimated the BH mass using the tight empirical cor- 
relation between BH mass and bulge stellar velocity dispersion 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000), adopting the 
Mbh-<^ relation given by Gebhardt et al. (2000) and a = 268 
km S-' from McElroy (1995)1 The VLBA flux densities of 7 
and 9 mJy at 1.7 and 5.0 GHz can be compared to the VLA 
peak flux densities of 13 and 10 mJy at similar frequencies in 
late 1999 (HUOl), so there may be weak 1.7-GHz emission on 
scales between tens of mas and 1 arcsec. 

NGC 4203.— NGC 4203 is a nearly face-on SO galaxy with 
a LINER 1.9 nucleus (Ho et al. 1997). It has a double-peaked 
broad Ha emission line with a full-width near zero intensity of 
at least 12,500 km s"' (Shields et al. 2000), and a fairly weak 
nuclear X-ray source detected with Chandra, having Lx(2- 
10 keV) = 5.0 X lO^** erg s"' (Ho et al. 2001). For a = 124 
km S-' (DaUe Ore et al. 1991), Mbh = 2.0 x 10^ Mq. The upper 
limit for the source radius is ~0.02 pc, or 1.0 x 10^/?s. and the 
VLBA flux densities are consistent with the 1999 VLA peaks 
(HUOl). 

NGC 4579.— NGC 4579 is an SABb galaxy with a type 1.9 
Seyfert or LINER nucleus (Ho et al. 1997). Its broad Ha hne is 
double-peaked or double- shouldered (Earth et al. 2001), and it 

contains a variable hard X-ray source (Terashima et al. 2000b) 
with a 2-10 keV luminosity of ~ 8.9 x lO'*" erg s"' as measured 
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TABLE 2 
Properties of the Radio Cores 



Galaxy 


D 


Mbh 






logP8.4 


log LsA 


log Lx 


log Lsa/Lx 


log Lbol/iEdd 




(Mpc) 


(Me) 






(W Hz-1) 


(erg s"^) 


(erg s"^) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 3147 


40.9 


3.6x10** 


0.20 ±0.10 


0.11 ±0.21 


21.28 


37.98 


41.53 


-3.55 


-4.30 


NGC 4203 


9.7 


2.0x10^ 


0.38 ±0.10 


0.09 ±0.21 


20.05 


36.76 


39.69 


-2.93 


-4.42 


NGC 4579 


16.8 


5.0x10'' 


0.20 ±0.10 


-0.31 ±0.21 


20.81 


37.52 


40.95 


-3.43 


-3.80 



Note. — Col. (1) Galaxy name. Col. (2) Distance from TuUy 1988, who assumes Ho = 75 km s^^ Mpc^^. Col. (3) BH mass estimated 
from the Mbh — c relation (see text). Col. (4) Two-point spectral index between 1.7 and 5.0 GHz. Col. (5) Two-point spectral index 
between 5.0 and 8.4 GHz. Col. (6) Monochromatic power at 8.4 GHz. Col. (7) Spectral luminosity at 8.4 GHz, Lsa = uPsa- Col. (8) 
X-ray (2-10 keV) luminosity (see text). Col. (9) Ratio of radio to X-ray luminosity. Col. (10) Ratio of bolometric luminosity (see text) 
to the Eddington luminosity. 



with Chandra (Ho et al. 2001). The BH mass inferred from the 
Mbh-<J relation is 5 x 10^ Mq (Barth et al. 2001). Our 
VLBA 8.4-GHz beam of 2.4 mas x 1.0 mas limits the source 
radius to ^-0.03 pc, or ~ 6200/?s- The 1.7-GHz VLBA flux den- 
sity is consistent with the previous VLA peak (HUOl), while 
the 5 -GHz value is considerably lower (22.8 mJy vs. 38 mJy). 
Falcke et al. (2001) found this galaxy to be variable at 15 GHz 
on time scales of 1-3 yr, so the change in the 5-GHz flux density 
over 17 months is most likely due to variabiUty. This variability 
implies that the core spectral index changes over time, so model 
interpretations based solely on the spectral index must be made 
with caution. 

4. IMPLICATIONS AND SUMMARY 

What is the emission mechanism of the radio cores? In 
the Seyfert galaxy NGC 1068, Galhmore, Baum, & O'Dea 
(1997) resolved the flat-spectrum component SI at 8.4 GHz us- 
ing the VLBA; they attributed S 1 to free-free emission from a 
parsec-scale torus. However, other flat-spectrum Seyfert cores 
are unresolved by the VLBA and appear to have other origins 
(Mundell et al. 2000). The free-free interpretation also is un- 
tenable for our sources because (1) they are much more com- 
pact (r ^ 0.05 pc) than the dimensions of nuclear tori and (2) 

they have brightness temperatures {Tb ^ lO' K; Table 1) more 
than two orders of magnitude above the values expected for 
thermal emission. The relatively flat spectra over a factor of 5 
in frequency (Fig. 2) seem to rule out free-free absorption as 
well, though NGC 4579 shows a sUght spectral turnover above 
5 GHz. Therefore, the prevalence of flat or inverted spectra in 
LLAGNs cannot, in general, be due to free-free emission or ab- 
sorption. The brightness temperature limits are consistent with 
classical synchrotron self-absorption, although more than one 
self-absorbed component would have to be present, since self- 
absorbed synchrotron radiation from a single power-law distri- 
bution of electrons should result in a « 2.5. 

If our galaxies host massive BHs, as suggested by current 
BH demography studies and the evidence for nonstellar activ- 
ity summarized in § 3, they must be radiating significantly be- 
low their Eddington limits. The nuclear bolometric luminosi- 
ties of NGC 4203 and NGC 4579 are known from Ho et al. 
(2000) and Ho (1999), respectively. For NGC 3147, we assume 
^-boi = 6.7Lx(2-10 keV), empirically derived from an average 
of 10 objects analyzed by Ho (1999) and Ho et al. (2000). As 
shown in Table 2, all three objects are highly sub-Eddington 

' The apparent spectral turnover in NGC 4579, however, would require some 
Narayan 1999). 



systems: Lboi/£-Edd < 10""*. In particular, they all lie comfort- 
ably within the expected value for ADAFs, which is estimated 
to be Lboi/iEdd < 10"^ (e.g., Narayan et al. 1998; Quataert 
2001). Additional arguments for the existence of ADAFs in 
NGC 4203 and NGC 4579 are given in Ho et al. (2000) and Ho 
(2001), among them being the absence of the optical/ultraviolet 
continuum bump normally attributed to thermal emission from 
a thin disk. Thus, the compact radio cores could be emission 
from ADAFs. This is consistent with (1) the compactness of 
the sources (r < lO'^.Rs), (2) the apparent self-absorbed syn- 
chrotron nature of the emission, (3) the brightness temperature 
limits, and (4) the spectral indices^. 

Notwithstanding these arguments, the majority of the radio 
emission most likely does not come from accretion flows. We 
arrive at this conclusion by comparing the luminosity output of 
the radio and the X-ray bands. Table 2 shows that the ratio of 
the spectral luminosity at 8.4 GHz to the luminosity in the 2-10 
keV band spans 3 x lO""* <,Ua/Lx <,! x 10"^ Yi & Boughn 
(1998) give a convenient expression for the relation between ra- 
dio and X-ray luminosities in a standard ADAF model. Using 
our values of BH mass and X-ray luminosity in their Equation 
2.10, we predict values for L^a/Lx that are ~10 times smaller 
than the observed values. Hence our sources are overluminous 
in the radio compared to the X-rays. We note that this discrep- 
ancy is not alleviated by appealing to more recent variants of 
the ADAF model that incorporate outflows or convection (e.g., 
Quataert 2001). Either of these effects will suppress the ra- 
dio component with respect to the X-rays (Quataert & Narayan 
1999; Ball, Narayan, & Quataert 2001), exacerbating the prob- 
lem. If the radio core were larger than in the canonical ADAF 
models, the surrounding X-ray bremsstrahlung region would 
be correspondingly larger, again increasing the discrepancy be- 
tween our high observed values of L^a/Lx and the lower model 
values. 

The most likely origins for the radio cores in our three galax- 
ies are compact jets or outflows, as discussed by Falcke (2001) 
and Nagar et al. (2001). The high brightness temperatures, the 
flat/inverted spectra, and the compactness of the cores easily 
can be reproduced in such a model. However, it is unclear if the 
jet model can self-consistently explain both the radio and X-ray 
emission, let alone the entire spectral energy distribution. Fal- 
cke & Markoff (2000) successfully fitted the radio and X-ray 
spectra of Sgr A* with the jet model, but L^a/Lx in Sgr A* 

modification of the standard ADAF model, such as inclusion of winds (Quataert & 
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Fig. 2. — VLBA radio spectra of the unresolved cores of the three LLAGNs. 
The spectrum of the Galactic Center source Sgr A* (Falcke et al. 1998) is 
shown for comparison, assuming a distance of 8.0 kpc (Reid 1993), with power 
scaled up by a factor of 1000. Error bars for the powers are smaller than or 
equal to the symbol size. 



is ^100 times higher than in our objects during the quiescent 
X-ray phase found by Baganoff et al. (2001). For a given radio 
jet power, the X-ray luminosities from self-Compton emission 
are well-constrained (e.g., Falcke & Markoff 2000), and are un- 
likely to contribute significantly in our galaxies. 

A natural way to account simultaneously for the full spectral 
energy distribution of LLAGNs, such as those discussed here, is 
to incorporate facets of both classes of models discussed above, 
namely a jet and an ADAF (or some closely related variant of 
a low radiative efficiency accretion flow; see Quataert 2001). 
This is the approach taken by Yuan (2000), but we note that no 
existing model self-consistently incorporates both components. 

We thank Heino Falcke for helpful discussions and Nadia 
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ules. This research has made use of NASA's Astrophysics 
Data System Abstract Service and the NASA/IPAC Extragalac- 
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Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. The 
work of L. C. H. is partly funded through NASA grants awarded 
by the Space Telescope Science Institute, which is operated by 
AURA, Inc., under NASA contract NAS5-26555. 
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